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EXECUTIVE  SUMMARY 


A  need  exists  to  be  able  to  detect  overflying  satellites  reliably  from  the  ground,  whether  for  space 
surveillance,  signature  studies,  or  antisatellite  systems.  Infrared  sensors  are  suggested  as  a  possible  means; 
they  are  passive  and  offer  much  belter  angular  position  determination  than  a  microwave  radar.  Funher- 
more,  an  infrared  system  is  much  smaller  and  consumes  less  power  than  a  radar.  Clear  weather  is  required 
for  successful  operation,  but  it  may  be  feasible  to  place  the  infrared  sensor  in  an  aircraft  to  fly  above 
cloud  cover. 

Optical  properties  of  the  atmosphere  impose  limitations  on  any  infrared  sensor  system  for  satellite 
acquisition.  The  important  effects  are  attenuation  of  the  satellite  radiation  signal,  infrared  radiance  from 
the  atmosphere  itself,  and  turbulence  which  blurs  the  image  of  the  satellite.  We  have  used  the  LOWTRAN7 
code  to  model  the  transmission  and  radiance.  Two  wavelength  bands,  3-5  pm  and  8-12  pm,  provide 
windows  where  the  transmission  is  high  enough  and  radiance  low  enough  for  direct  observation  of 
satellites.  Image  spreading  by  turbulence  decreases,  but  spreading  by  diffraction  increases  with  increasing 
wavelength.  For  a  modest-sized  (24-inch)  telescope,  the  combined  spreading  effect  is  at  its  minimum 
around  5  pm  in  wavelength,  depending  on  the  elevation  angle. 

The  range  a:  which  a  given  satellite  may  be  detected  varies  as  the  founh  root  of  the  product  of  the 
detector  quantum  efficiency  and  the  number  of  detecting  elements  in  the  array.  Only  recently  have 
infrared  focal-plane  arrays  become  available  which  have  enough  delecting  elements  and  detection  effi¬ 
ciency  to  be  effective  in  acquiring  most  satellites.  Large  detector  arrays  can  now  be  obtained  for  both  the 
3-5  pm  mid-wave  and  the  8-12  pm  long-wave  infrared  bands.  We  have  evaluated  the  potential  perfor¬ 
mance  of  actual  arrays  based  on  InSb.  HgCdTc,  Si.Ga,  PiSi,  and  IrSi  detectors. 

The  performance  calculations  of  this  report  assume  that  the  infrared  sensor  array  is  applied  to  a 
reflecting  telescope  of  60-cm  aperture.  Uncooled  optics  are  used  to  relay  the  satellite  image  to  the  focal 
plane.  The  relay  optics  also  serve  to  re-image  the  primary  mirror  at  the  cold  apenure  stop  in  the  detector 
dewar.  Background  radiation  from  the  telescope  and  its  surroundings  is  thus  minimized.  We  estimate  that 
an  optical  transmission  of  76  percent  may  be  achieved.  A  cold  filter  in  front  of  the  focal-plane  array 
passes  only  those  wavelengths  for  which  the  atmosphere  is  most  transparent.  The  optical  system  must 
include  a  means  to  calibrate  the  responsivities  and  offsets  of  the  individual  detector  elements.  A  flat-field 
illuminator  is  easily  provided.  Detector  data  are  recorded  at  more  than  one  radiance  level,  bracketing  the 
level  expected  during  an  acquisition  search.  A  set  of  coefficients  derived  from  the  calibration  data  are 
applied  to  correct  the  real-time  detector  data. 

Immediately  prior  to  an  acquisitic  i  search,  atmospheric  background  data  are  recorded  digitally.  This 
will  require  considerable  precision,  because  the  background  radiance  is  very  large  relative  to  the  satellite 
signals  to  be  detected.  The  procedure  for  acquisition  is  to  choose  an  elevation  angle  at  which  the  satellite 
can  be  detected  with  confidence,  and  scan  repeatedly  across  the  nominal  trajectory  of  the  satellite  until 
it  is  detected.  Each  scan  is  likely  to  consist  of  several  integration  periods,  during  which  the  telescope 
tracks  the  nominal  angular  peed  of  the  satellite,  separated  by  rapid  steps  to  the  next  field  of  view.  During 
the  “staring"  period,  the  frames  from  the  infrared  camera  are  summed  to  provide  the  maximum  feasible 
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integration  time.  This  is  essential  to  achieve  the  best  possible  signal-to-noise  ratio  (SNR).  If,  after 
subtracting  the  corresponding  background  signal,  a  spot  image  remains  with  a  net  signal  over  six  limes 
the  noise  level,  an  object  detection  is  declared  and  the  telescope  begins  tracking  the  candidate  object. 

For  each  of  the  candidate  sense r  arrays  the  satellite  acquisition  performance  has  been  calculated, 
assuming  for  simplicity  that  'he  satellite  is  in  a  circular  orbit  passing  overhead.  It  is  assumed  that  the 
nominal  satellite  orbital  elemei'ts  arc  known,  but  the  present  position  is  uncenain  to  within  a  radius  of 
50  km,  perhaps  ab  a  result  ot  a  recent  maneuver.  The  allowable  integration  times  in  such  a  case  depend 
on  where  the  satellite  is  at  the  time  of  observation.  For  the  typical  system  parameters  used  here,  they  are 
of  order  0.1-1  second.  Calculations  of  transmitted  signals  and  radiant  backgrounds  assumed  midlatitude 
summer  conditions.  Rural  haze,  23-km  visibiliiy,  and  observations  from  sea  level  are  assumed.  The 
dominant  noise  is  from  the  background  radiances,  both  atmospheric  and  instrumental.  However,  reported 
values  of  minimum-irradiance  system  electronic  noise  have  been  included  in  the  analysis. 

The  results  of  the  calculations  may  be  seen  in  two  ways.  In  the  first,  a  nominal  satellite  size  and 
temperature  are  chosen,  and  the  zone  of  range  and  elevation  within  which  each  type  of  infrared  camera 
can  detect  the  satellite  is  plotted.  We  find  that  the  HgCdl  e  vsiem  can  detect  a  5-m’,  3(X)-K  satellite  almost 
1 1,000  km  away.  The  Si:Ga  system  can  detect  it  a^  km.  Of  the  sensors  operating  in  the  3-5  pm 
band,  InSb  performs  best,  with  over  5,000-km  detection  range.  IrSi  and  PtSi  systems  can  detect  at  2,000 
and  1,200  km,  respectively.  These  ranges  bold  for  elevations  above  45®.  Closer  to  the  horizon,  increasing 
atmospheric  absorption  and  radiance  cause  the  detection  range  to  fall  off  rapidly. 

The  second  type  of  presentation  is  a  set  of  contour  plots  of  minimum  deicciabie  infrared  signature 
as  a  function  of  satellite  position.  On  the  plot  for  each  detector,  a  grid  is  superimposed  to  show  circular 
orbits  at  various  altimdes  and  lines  of  fixed  elevation  angle.  Given  the  altitude  and  signature  of  a  satellite, 
the  range  and  elevation  at  which  it  could  be  detected  by  each  camera  is  readily  found. 

We  conclude  that  all  the  proposed  systems  may  be  used  to  acquire  satellites  in  low  earth  orbits. 
Whether  a  given  detector  system  will  fulfill  a  particular  requirement  will  depend  on  the  actual  altitudes 
and  infrared  signatures  of  the  satellites  sought.  While  the  present  results  are  very  encouraging,  we 
consider  it  highly  desirable  to  validate  these  calculations  with  an  actual  experimental  test.  The  calcula¬ 
tions  arc  based  on  considerable  signal  averaging  and  background  subtraction.  Noise  effects  not  obeying 
simple  Gaussian  statistics,  such  as  turbulent  sky  radiance  fluctuations  or  "l/f’  noise  in  the  detectors, 
could  potentially  set  more  stringent  limits  to  the  perfonmance  of  an  infrared  satellite  acquisition  system. 
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1.  INTRODUCTION 


Ll  GROUND-BASED  DETECTION  OF  SATELLITES 

The  capability  to  detect  overflying  satellites  from  the  ground  has  a  number  of  applications,  including 
space  surveillance,  signature  studies,  and  antisatellite  systems.  The  means  to  acquire  the  satellites  include 
radar  and  passive  electro-optical  (E-0)  sensors.  Microwave  radar  systems  are  well  developed,  but  to  find 
a  typical  satellite  at  a  range  of  several  thousand  kilometers  they  must  employ  very-high-power  transmit¬ 
ters.  We  have  examined  the  possibility  of  acquiring  satellites  from  the  ground  using  the  infrared  sensors 
currently  available.  The  potential  performance  of  Schottky-barrier  infrared  sensors  for  the  surveillance 
application  has  been  studied  by  Cantella  [1].  That  report  primarily  considered  open  searches  for  previ¬ 
ously  unknown  objects  in  space.  It  is  assumed  in  this  study  that  the  elements  of  the  satellite  orbit  are 
known  approximately,  but  that  the  position  of  the  satellite  is  uncertain  to  within  a  specified  radius.  This 
report  examines  the  potential  performance  of  several  different  kinds  of  infrared  sensors  in  the  satellite 
acquisition  role. 

1.2  ADVANTAGES  OF  INFRARED  DETECTION 

The  infrared  (IR)  sensors  we  have  in  mind  would  be  used  to  passively  detect  the  radiant  emissions 
of  the  satellite.  Unlike  radar,  a  passive  sensor  system  does  not  advertise  its  presence  and  location.  This 
would  be  a  great  advantage  for  a  mobile  military  system;  to  protect  a  fleet  from  observation,  for  instance. 
An  IR  sensor  requires  a  telescope  of  modest  size  (~24-inch  aperture)  to  gather  radiation  from  the  satellite; 
the  same  sensor  can  provide  necessary  tracking  information  after  acquisition.  The  angular  precision 
attained  by  the  IR  acquisition  system  should  be  much  better,  of  order  1 0  prad,  as  opposed  to  1  mrad  for 
a  radar  acquisition  system.  A  radar  system  for  satellite  acquisition,  with  its  associated  power  supplies,  is 
large  and  heavy.  In  contrast,  it  may  be  feasible  to  install  an  infrared  acquisition  system  in  an  aircraft, 
which  can  fly  above  the  weather  (this  may  be  required  in  some  applications,  as  the  IR  system  requires 
clear  sky  to  operate). 

Given  the  advantages  of  passive  E-O  sensors,  the  next  question  is  wl;ich  sensor  is  best  Visible 
detectors  are  simple  to  use  and  work  very  well  in  the  terminator  mode  (satellite  in  sunlight,  observer  in 
darkness).  However,  it  is  desirable  to  be  able  to  acquire  a  satellite  at  any  time,  day  or  night.  To  find  a 
satellite  in  the  earth's  shadow,  it  is  necessary  to  sense  its  infrared  emissions.  There  are  two  wavelength 
ranges  in  the  IR  in  which  the  atmosphere  is  relatively  dark  and  transparent,  even  in  the  daytime.  They 
are  called  midwave  infrared  (MWIR)  (3-5  pm),  and  long- wave  infrared  (LWIR)  (8-12  pm). 

13  OBSTACLES  TO  INFRARED  ACQUISITION 

Until  very  recently,  imaging  IR  sensors  (sometimes  called  FLIRs)  relied  on  one  or  a  small  number 
of  detectors  optically  scanned  over  a  scene  to  produce  an  image.  Such  devices  arc  unlikely  to  be  sensitive 
enough  to  detect  a  satellite  at  the  required  ranges  because  the  detector  dwells  for  only  a  short  time  on 
any  one  picture  element  (pixel).  However,  several  types  of  staring  focal-plane  array  (FPA)  sensors  are 
now  available  for  the  IR  bands.  They  include  Schottky-barrier  photovoltaic  arrays  on  silicon  CCDs  and 
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hybrid  arrays  of  InSb,  Hgi-^Cd^T-e,  or  Ga-doped  silicon  on  silicon  complementary  metal-oxide  semi¬ 
conductor  (CMOS)  multiplexors.  The  large  number  of  detectors  in  these  staring  arrays  make  it  possible 
to  search  a  reasonable  solid  angle  with  the  necessary  sensitivity. 

The  atmosphere  not  only  absorbs  IR  radiation  to  some  extent  at  all  wavelengths,  it  also  emits  IR 
radiation  to  a  corresponding  degree.  This  is  true  even  in  the  MWIR  and  LWIR  “windows”  mentioned 
earlier.  For  some  proposed  applications,  it  is  necessary  to  acquire  the  satellite  at  the  lowest  possible 
elevation  angle,  where  the  atmospheric  transmission  and  radiance  effects  are  at  their  worst.  We  have 
attempted  to  analyze  the  potential  performance  of  a  practical  IR  acquisition  system,  using  currently 
available  sensors,  in  relation  to  actual  satellite  signatures  and  kno^  n  optical  properties  of  the  atmosphere. 

1.4  PROPOSED  EXPERIMENT 

An  experimental  test  of  IR  acquisition  of  satellites  would  be  highly  desirable  to  verify  our  perfor¬ 
mance  calculations,  make  further  measurements  of  satellite  signatures,  and  explore  the  practical  limita¬ 
tions  of  the  present  technology.  The  MIT  Lincoln  Laboratory  Firepend  telescope  is  an  ideal  testbed  for 
this  purpose,  as  it  is  already  well-equipped  to  track  satellites.  It  includes  an  auxiliary  telescope  of  60-cm 
aperture,  which  appears  suitable  for  this  purpose.  Therefore,  the  performance  calculations  in  this  report 
arc  based  on  the  use  of  that  telescope.  Comparison  with  other  apertures  is  very  straightforward. 

1.5  ORGANIZATION  OF  REPORT 

Sections  2  and  3  of  this  report  outline  the  characteristics  of  likely  satellites  and  the  relevant  properties 
of  the  atmosphere.  Section  4  describes  an  IR  acquisition  system  and  several  different  types  of  detectors 
which  it  might  use.  Section  5  presents  an  analysis  of  the  performance  which  may  be  expected  from  these 
detectors  based  on  fundamental  physical  limitations.  For  the  IR  acquisition  system,  the  key  result  is  the 
maximum  range  of  detection  for  a  given  IR  signature  as  a  function  of  satellite  elevation  angle. 
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2.  SATELLITE  CHARACTERISTICS 


2.1  INFRARED  SIGNATURES 

Aji>  object  emits  electromagnetic  radiation  by  virtue  of  its  thermal  energy.  The  emittance  is  given  by 
E  =  eoT"*  ,  (I) 

where 

e  =  emissivity  of  the  surface* 

(T  =  the  Stefan-Boltzmann  constant  =  5.67  x  10  ®  (W  m**  K"*) 

T  =  the  absolute  temperature,  in  Kelvins  (K). 

The  blackbody  radiance  within  the  MWIR  (3-5  pm)  and  LWIR  (8-12  pm)  wavelength  bands  is 
plotted  as  a  function  of  the  temperature  of  the  object  addressed  in  Figure  1 .  Figure  2  is  a  plot  of  the  ratio 
of  ihe  radiance  in  the  MWIR  band  to  that  in  the  LWIR  band  as  a  function  of  emitter  temperature.  Such 
a  plot  can  be  used  to  estimate  a  sati^llite  temperature  given  the  measured  signatures  in  the  two  bands. 


Figure  I.  Blackbody  radiance  within  the  MWIR  (3-5  pm)  and  LWIR  (8-12  pm)  spectral  bands  vs.  source 
temperature. 
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Figure  2.  Ratio  of  blackbody  radiance  within  the  MWIR  (3-5  \im)  band  to  that  in  the  LWIR  (8-12  pm)  band  vs. 
source  temperature. 

In  addition  to  the  thermal  radiation  emitted  by  a  satellite,  there  is  some  reflected  radiation.  Reflected 
sunlight  will  especially  increase  the  infrared  signature  in  the  3-5  pm  band.  The  tarth  will  bathe  the 
satellite  in  thermal  radiation  with  an  effective  temperature  approximately  equal  to  the  local  ground-level 
temperature.  A  fraction  of  this  earth  radiation  (1-e,  where  e  is  the  satellite  emissivity)  will  be  reflected 
back  down  toward  the  observer. 

For  the  satellites  of  interest  in  this  study,  we  assume  that  the  projected-area  enittance  products  range 
from  1-25  m^.  llie  satellite  surface  temperatures  should  be  roughly  280-300  K;  then  the  total  LWIR 
(8-12  pm)  signatures  are  in  the  range  25-1,000  W/ster.  Assuming  that  the  emittance  is  constant  over  the 
3-12  pm  wavelength  range,  the  corresponding  MWIR  (3-5  pm)  signatures  are  roughly  1-50  W/ster. 

2.2  TRAJECTORIES 

The  satellites  to  be  acquired  are  assumed  to  be  in  relatively  low  earth  orbits,  (250-2,500-km 
altitude).  The  nominal  orbital  elements  of  the  satellite  are  assumed  to  be  known;  however,  a  robust 
acquisition  system  should  be  able  to  acquire  a  satellite  which  has  just  maneuvered  to  alter  its  orbit.  The 


4 


(EARTH) 


SEARCH  BOX  ON  SKY 


da/dt 


NOMINAL 
ACQUISITION 
ELEVATION.  O 


POSITION 

UNCERTAINTY 

RADIUS 


Figure  3.  Geometry  of  acquisition  scenario  (sc  hematic  1 


resulting  position  uncertainty  zone  may  typically  be  50  km  in  radius.  To  simplify  the  analysis,  it  was 
assumed  that  the  nominal  satellite  orbit  is  circular  and  passes  directly  over  the  observer.  The  relative 
positions  are  diagrammed  in  Figure  3. 

For  some  systems,  llie  satellite  must  be  acquired  at  the  lowest  feasible  elevation  angle.  This  angle 
is  chosen  to  ensure  a  high  probability  of  acquisition  as  the  satellite  reaches  that  elevation.  The  acquisition 
telescope  is  scanned  across  the  position  uncertainty  zone  just  rapidly  enough  to  complete  the  search  as 
the  zone  rises  through  the  chosen  elevation.  The  orbital  velocities  of  satellites  in  low  earth  orbits  are 
roughly  7  km/sec.  Their  angular  speeds  as  seen  from  the  ground  will  range  between  1  mrad/sec  for  a 
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saiellile  just  rising  above  the  horizon  to  30  mrad/sec  for  a  very  low  satellite  passing  overhead.  If  acqui¬ 
sition  is  to  take  place  at  elevation  angles  less  than  30°  the  angular  rates  will  not  exceed  5  mrad/sec. 
Detailed  expressions  for  the  angular  speeds  of  satellites  are  presented  in  Appendix  B. 
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3.  ATMOSPHERIC  OPTICS 


3.1  RADIANCE  AND  TRANSMISSION  VS.  ELEVATION,  3-14  yim 

The  opacity  of  the  atmosphere  in  the  infrared  region  is  a  complicated  function  of  wavelength.  It 
depends  on  the  relative  concentration  of  water  vapor  and  aerosols,  which  vary  with  time  and  place. 
Fonunately,  the  problen.  has  been  studied  extensively,  and  a  computer  code,  LOWTRAN,  has  been 
developed  which  provides  a  numerical  model  of  the  atmosphere  suitable  for  our  analysis  [2].  We  have 
obtained  the  ground-to-space  IR  transmission  and  radiance  of  the  atmosphere  as  a  function  of  wavelength 
and  elevation  angle  by  running  LOWTRAN7  on  a  VAX  computer.  Conditions  assumed  for  these  runs 
were  middle  latitude,  summer,  and  rural  haze  aerosols  resulting  in  23-km  visibility.  These  conditions  are 
considered  representative  of  the  least-favorable  season  for  infrared  observations  of  space.  (An  IR  sensor 
will  perfonn  better  in  winter,  because  the  lower  winter  air  temperature  and  water  vapor  content  lead  to 
lower  atmospheric  opacity.)  The  transmission  is  plotted  in  Figure  4;  the  radiance  is  shown  in  Figure  5. 
Notice  that  the  wavelengths  of  greatest  transmittance  are  also  those  at  which  the  radiance  is  least,  as 
absorplance  and  radiance  are  directly  related.  At  the  wavelengths  of  least  transmittance,  the  radiance 
resembles  that  of  a  blackbody  at  295  K.  Scattered  solar  radiation  is  not  included  in  the  plotted  radiance, 
so  the  curves  are  most  appropriate  to  nighttime  observation.  This  radiation  equals  the  low-elevation 
thermal  radiance  at  approximately  3  pm  and  decreases  with  increasing  wavelength.  In  daytime,  the 
scattered  solar  radiation  will  increase  the  MWIR  background  radiation  by  a  fraction,  depending  on  the 
relative  angle  of  the  sun  from  the  look  direction. 


3.2  TURBULENCE 


Turbulence  in  the  atmosphere  can  also  affect  observations,  increasing  the  size  of  the  image  of  a  point 
object  which  would  otherwise  be  limited  by  diffraction.  (At  the  ranges  of  detection  anticipated,  the 
satellites  would  be  unresolvable.)  We  assume  that  the  optics  are  good  enough  at  these  wavelengths  to  be 
diffraction-limited.  The  angular  size  of  the  image  0(full  width  at  half-maximum,  m  radians)  is  then  given 
by 


02  = 


/  \2 
1.22A 

<  ^0  ; 


(2) 


where  A  is  the  wavelength,  D  the  diameter  of  the  telescope  objective;  the  atmospheric  coherence  length 
is  given  by  [3] 


rQ=0.185| 


A2 
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where  is  the  structure  constant  of  the  refractive  index  fluctuations. 


(3) 
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TRANSMISSION 


WAVELENGTH  (^m) 


Figure  4.  Atmospheric  transmission  vs.  wavelength  for  various  indicated  elevation  angles  ( in  degrees).  Midlatitude 
summer  conditions  with  rural  aerosol  haze  giving  23 -km  visibility  have  been  assumed. 
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Figure  5.  Atmospheric  radiance  vs.  wavelength  for  various  indicated  elevation  angles  (in  degrees).  Conditions 
assumed  are  the  same  as  for  Figure  4.  No  scattered  solar  radiation  is  included. 

In  Equation  (2),  the  first  term  is  the  contribution  of  diffraction  to  the  image  size;  the  second  is  caused 
by  turbulence.  If  we  know  the  coherence  length  or,  alternatively,  the  size  of  the  atmospheric  “seeing  disk” 
for  visible  light,  we  need  only  the  wavelength  dependence  given  in  Equation  (3)  to  find  the  value  of  rQ 
for  infrared  wavelengths.  It  is  then  unnecessary  to  know  the  structure  constant  for  the  whole  atmospheric 
path.  It  can  also  be  shown  from  Equation  (3)  that  the  dependence  of  tq  on  the  elevation  angle,  a,  is 

r^  =  r^|  (sina)^^^  .  (4) 

0  0lct=90“ 

Typical  values  of  roat  high  elevations  for  visible  light  {X  *  0.5  pm)  arc  5  cm  in  daytime  conditions, 
increasing  to  10  cm  in  good  nighttime  conditions.  If  we  call  this  value  rQg£F  ,  then  we  may  re -express 
Equation  (2)  as 


where  Xo  is  the  wavelength  for  rQg£f  . 
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For  ihe  system  under  analysis  we  have  O  =  60  cm;  we  will  assume  rc^EF  to  be  5  cm  {-Ao  =  0.5  |im) 
at  a  =  90°.  With  these  assumptions,  the  separate  contributions  of  diffraction  and  turbulent  seeing  con¬ 
ditions  to  the  size  of  an  image  are  plotted  against  wavelength  (for  various  elevation  angles)  in  Figure  6. 
The  resulting  total  image  size  is  plotted  in  Figure  7. 

We  find  that  for  the  low  elevation  angles  atmospheric  turbulence  is  a  very  significant  contributor  to 
the  size  of  the  satellite  images.  In  the  MWIR  (3-5  pm)  band  it  is  the  main  contributor,  while  for  the 
LWIR  (8-12  pm)  band  it  is  at  least  comparable  to  diffraction. 


Figure  6.  Contributions  of  diffraction  and  turbulence  ("seeing")  to  image  size  vs.  wavelength.  Turbulence  cunes 
are  for  elevation  angles  5-90°.  Diffraction  size  is  for  dO-cm  (24 -inch)  objective  diameter. 
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IMAGE  SIZE  (prad) 


Figure  7.  Total  size  of  point-source  image  vs.  wavelength  for  elevation  angles 
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4.  SYSTEM  DESKJN 


The  esscniial  elenienis  ol  a  system  for  infrared  aeqiiisiiion  and  tracking  of  satellites  are  shown  in 
Figure  H.  Detads  ol  the  special  subs\ stems  and  operating  procedures  are  discussed  in  the  tollouing 
sec  lions. 


hijfurc  S.  Schematic  diagram  of  an  infrared  satellite  aiywsition  system 


4.1  OPTICAL  ELEMENTS 

The  numerical  calculations  which  follow  are  based  on  use  of  a  60-cm  Newtonian  telescope  such  as 
the  auxiliary  telescope  of  the  Firepond  Obsenatory.  As  nia>  be  seen  from  Figure  6.  this  aperture  is 
sufficient  to  give  diffraction-limited  images  comparable  to  or  smaller  than  the  atmospheric  seeing  disk 
Additional  optics  will  be  included  to  provide  a  suitable  effective  focal  length  for  the  focal  plane  array 
and  to  ensure  effective  cold-stopping  and  baffling  to  reject  stray  radiation. 
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4.1.1  Primary  and  Diagonal  Mirrors 

The  primary  mirror  has  a  diameter  of  0.60  meters  and  a  focal  length  of  3.25  meters.  The  collecting 
area  is  0.27  m“  after  allowance  for  obscuration  by  the  diagonal  flat  mirror. 

4.1.2  Relay  Lenses 

A  relay  lens  assembly  provides  the  r  esired  effective  focal  length.  The  total  field  of  view-  of  the  FPA 
should  be  large,  to  allow  maximum  inir ^ration  time  for  each  possible  satellite  position;  this  suggests  the 
use  of  the  shortest  practical  effective  focal  length.  Tne  effective  focal  length  is  cho.sen  so  that  the  image 
size  is  approximately  one  pixel.  Detectors  much  larger  than  the  satellite  image  will  receive  more  back¬ 
ground  radiance  than  is  necessary,  lei  ding  to  a  decrease  in  the  SNR.  Furthermore,  if  the  .satellite  image 
were  much  smaller  it  could  sometimes  fall  on  the  nonresponsive  spaces  between  detectors. 

A  longer  focal  length  would  lead  to  greater  precision  in  tracking  after  initial  acquisition.  The  SNR 
need  not  be  adversely  affected  by  an  increase  in  image  size  to  several  pixels,  provided  that  the  system 
electronic  noise  and  readout  noise  are  small  compared  to  the  background  photon  noise,  i.e.,  the  detector 
is  background-limited.  However,  because  the  present  study  is  primarily  concerned  with  acquisition,  we 
will  assume  that  maximum  field  of  view  is  the  most  important  consideration. 

For  operation  at  low  elevation  angles,  the  image  size  typically  will  be  25  |irad  for  all  IR  wavelengths 
of  interest  (see  Figure  7).  The  typical  pixel  pitch  of  the  available  FPAs  is  50  |im.  Thus,  the  effective  focal 
length  should  be  about  2  meters,  which  gives  an  effective  f/number  for  the  60-cm  telescope  of  f/3.35.  The 
magnification  of  the  relay  optics  is  2/3.25,  or  0.615. 

4.U  Cold  Stop  and  Baffling 

The  introduction  of  the  relay  lens  provides  a  simple  way  to  locate  a  cold  stop  so  that  only  desired 
radiation  from  the  optical  aperture  is  received  at  the  focal  plane.  An  image  of  the  primary  mirror  appears 
between  the  relay  lens  and  the  focal  plane.  This  is  to  be  located  inside  the  dewar;  a  cold  stop  is  placed 
there.  The  cold  stop  structure  is  to  include  a  central  obscuring  disk  located  at  the  image  of  the  diagonal 
mirror.  Other  cold  masking  elements  block  radiation  from  the  diagonal  supports. 

Cold  baffling  elements  in  the  dewar  should  be  designed  to  trap  stray  radiation  before  it  can  be 
scattered  onto  the  FPA.  Objects  outside  the  dewar  which  are  above  ambient  temperature  should  be 
shrouded  from  the  optical  path. 

4.1.4  Cold  Filter 

Depending  on  detector  characteristics,  a  cold  passband  filter  should  be  included  in  the  dewar  to  tailor 
the  system  response  to  just  those  wavelengths  which  are  transmitted,  not  radiated,  by  the  atmosphere. 
Simulation  calculations  have  been  performed  to  choose  the  best  passband  to  maximize  the  SNR. 
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4.1.5  Optical  Transmission 

We  will  assume  the  reflectivity  of  the  telescope  mirrors  to  be  95  percent  and  emissivity  to  be  5  percent. 
Refractive  optical  elements,  such  as  germanium  lenses  with  multilayer  anti  reflect  ion  coatings,  are  expected 
to  transmit  95  percent.  If  we  allow  for  three  elements  in  the  relay  lens,  plus  the  dewar  window,  we 
estimate  that  the  overall  optical  transmission  Xq  =  74  percent.  (The  cold  filter  is  treated  separately.) 

4.1.6  Emission  of  Optical  Elements 

Appendix  A  shows  that  the  effective  thermal  emissivity  of  the  optical  elements  at  ambient  tempera¬ 
ture.  as  .seen  by  the  FPA,  is  given  approximately  by 

n 

^SFF=  '  -  n^.  =  '  - 
1 

Thus,  given  the  assumptions  above,  we  find  the  value  of  0.26  for  e£/r/r  At  295  K,  the  total  blackbody 
radiance  from  the  optical  surfaces  is  thus  36  Wm'^ster  *. 

4.1.7  Scattered  Light 

Another  source  of  background  radiation  is  ambient  infrared  light  scattered  from  the  optical  surfaces. 
We  would  need  to  know  the  bidirectional  reflectance  distribution  function  (BRDF)  of  each  surface  to 
evaluate  this  fully.  For  this  study,  we  will  assume  that  the  condition  of  the  surfaces  is  good  enough  that 
this  background  is  small  compared  to  the  self-emission  from  the  surfaces. 

4.2  DETECTORS 

Characteristics  of  a  variety  of  IR  FPA  detectors  which  are  commercially  available  or  soon  to  be 
available  are  listed  in  Table  1.  The  indium-antimonide  and  mercury-cadmium-telluride  detectors  are 
arrays  of  photodiodes  connected  to  readout  multiplexor  arrays  by  indium-bump  bonds.  They  are  available 
from  Amber  Engineering,  Inc.  and  other  manufacturers.  The  parameters  of  these  detector  arrays  are  based 
on  those  of  a  128  x  128-clement  array  already  available  from  Amber  Engineering  (4).  Also  developed  by 
Amber  is  the  gallium-doped  silicon  array  with  a  similar  readout,  which  operates  as  a  photoconductor.  For 
comparison,  we  list  platinum-silicide  and  iridium-silicide  Schottky- barrier  detectors.  Although  they  have 
relatively  low  quantum  efficiencies,  their  monolithic  silicon  construction  and  CCD  readout  structures 
offer  high  uniformity,  good  pixel  yield,  and  low  noise.  The  properties  listed  for  the  PtSi  Schottky  array 
are  based  on  those  of  PtSi  cameras  produced  by  David  Samoff  Research  Center  [5].  The  data  for  the  IrSi 
array  are  based  on  detectors  produced  jointly  by  Lincoln  Laboratory  and  Ford  Aerospace  [6,7,8].  Ap¬ 
proximate  responsivities  of  the  detectors  listed  in  Table  1  are  plotted  in  Figure  9.  The  mociel  parameters 
have  been  adjusted,  where  possible,  to  match  measured  characteristics. 
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TABLE  1 


Properties  of  Infrared  Focal  Plane  Arrays 


Type  of  detector 

InSb 

HgCdTe 

Sl:Ga 

PtSl 

IrSI 

Mode 

PV 

PV 

PC 

PV 

PV 

Cutoff  wavelength  (pm) 

5.5 

10.5 

19.1 

~  5 

94 

Operating  temp.  (K) 

80 

55 

20 

80 

40 

Rows  X  columns 

256  X  256 

128  X  128 

128  X  128 

244  X  160 

256  X  256 

Pixel  pitch  (pm) 

50 

50 

50 

80 

40 

Filling  factor 

~  90  % 

>  64  % 

~  100  % 

39  % 

52  % 

Quantum  efficiency 

53  % 

-  40  % 

-  12  % 

~  1  % 

~  1  % 

Readout  noise  (e  /pixel) 

800 

1,600 

1,600 

300 

300 

Max.  charge  (eVpixel) 

4x10^ 

4x10^ 

4x10^ 

1.2x10® 

1.2x10® 

Good  pixel  fraction 

99  % 

<  97  % 

>  97  % 

99.9  % 

99.9  % 

Uniformity 

3.5  % 

-  5% 

-  5% 

-  2% 

-  2% 

4J  CALIBRATION  SOURCES 

The  responsiviiics  of  the  detectors  in  a  staring  FPA  typically  vary  by  about  3-*‘ ;  ctcem.  In  addition, 
there  are  various  sources  of  DC  offset  from  pixel  to  pixel.  It  is  essential  to  correct  the  output  of  the  array 
for  these  gain  and  offset  differences.  For  this  purpose,  a  flat-field  calibration  source  (or  sources)  should 
be  included  in  the  camera  optics  assembly.  For  calibration,  a  fold  mirror  in  the  optical  train  is  rotated  to 
a  position  such  that  radiation  from  the  calibration  source  illuminates  the  focal  plane  uniformly.  Just  prior 
to  each  observation,  camera  signals  are  to  be  recorded  for  at  least  two  different  levels  of  illumination. 
These  levels  should  closely  bracket  the  acmal  illumination  received  from  the  sky  during  an  acquisition 
search.  An  offset  and  gain  factor  are  computed  for  each  pixel;  these  arc  then  applied  to  the  raw  digital 
readouts  in  real  time.  This  is  called  a  “two-point”  correction  if  two  illumination  levels  are  used. 

4.4  CAMERA  INTERFACE  ELECTRONICS 

The  stream  of  data  from  any  of  the  IR  FPAs  considered  in  this  smdy  is  voluminous,  partly  because 
each  array  contains  lO^-lO^  pixels  in  each  frame.  Furthermore,  the  relatively  high  sky-background  radi¬ 
ance  and  limited  charge-holding  capacity  of  each  pixel  require  the  frame  lates  to  be  high  (30-1,000 
frames  per  second).  If  each  pixel  is  digitized  to  12  bits,  the  gross  data  rate  is  -4  x  10^  to  -10^  bits  per 
second. 
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Figure  9.  Numerically  modeled  responsivities  vs.  wavelength  for  some  candidate  FPA  detectors. 


For  the  acquisition  of  known  salelliies  the  telescope  can  track  the  nominal  positions  of  the  satellites. 
Thus,  the  satellites  are  nearly  stationary  within  the  field  of  view  and  a  long  integration  time  is  permissible 
before  the  image  data  are  sent  to  a  processor  for  satellite  detection.  The  sensitivity  of  the  search  increases 
with  integration  time. 

From  this  information  it  appears  advantageous  to  include  a  front-end  processor  in  the  camera  inter¬ 
face  electronics  to  integrate  the  image  data  as  it  comes  from  the  focal  plane.  At  its  minimum,  the 
processor  would  receive  digital  data  at  the  camera  frame  rate,  sum  the  frames  at  a  time,  and  produce 
output  data  at  \/N  of  the  camera  frame  rate.  Other  functions  should  be  included,  such  as  pixel-by-pixel 
correction  for  responsivity  and  offset  variations  and  provision  of  a  standard  composite  video  signal  for 
real-time  displays.  The  data  summation  provided  by  the  front-end  processor  must  not  be  confused  with 
data  sampling  to  provide  a  lower  frame  rate.  The  latter  process  does  not  provide  the  long  integration  times 
which  are  needed  here.  An  output  rate  of  30  frames  per  second  or  less  should  be  manageable  in  the 
subsequent  processors.  The  interface  electronics  should  also  include  scan-conversion  circuitry  to  produce 
a  standard  composite  video  signal  for  real-time  displays.  In  addition,  commercial  video-processing 
equipment  may  be  used  to  advantage. 
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4.5  HOST  COMPUTER 


The  acquisition  system  must  include  a  host  computer  (or  computers)  to  operate  effectively  in  real 
time.  Tasks  handled  by  the  host  include  pointing  the  telescope  toward  the  satellite,  controlling  its  tracking 
rates,  and  processing  the  IR  images  to  sense  the  presence  of  the  satellite  above  the  background  noise. 
Other  tasks  may  be  performed,  such  as  automatic  camera  calibrations  and  interactions  with  the  system 
operators.  However,  real-time  image  displays  should  also  be  provided  so  that  human  operators  may 
monitor  the  operations.  An  alert  human  can  probably  make  a  better  judgment  of  the  reality  of  a  suspected 
satellite  detection  than  a  computer  algorithm  can. 

4.6  ACQUISITION  PROCEDURE 

For  this  study  it  is  assumed  that  the  orbital  parameters  of  the  satellite  are  approximately  known.  The 
satellite  will  be  within  a  given  radius  of  the  predicted  nominal  position  (see  Figure  3).  This  uncertainty 
radius,  and  the  range  vector  to  the  nominal  position,  defines  an  angular  search  box  which  moves  with 
the  satellite  and  grows  as  the  satellite  approaches.  On  the  basis  of  the  expected  IR  signature  of  the 
satellite,  an  elevation  angle  is  chosen  which  gives  a  high  probability  of  detection  in  the  initial  search. 

Sky -background  data  are  recorded  immediately  prior  to  an  acquisition  observation.  This  would  be 
done  at  the  same  elevation  and  azimuth  angles  as  the  nominal  acquisition  search.  Thus,  the  level  and 
spectrum  of  sky  radiance  should  be  the  same  as  for  the  actual  search.  The  averaged  background  data  are 
then  subtracted  from  the  averaged  search  frames. 

The  box  is  to  be  searched  at  a  nearly  constant  elevation  angle  by  scanning  the  camera  in  azimuth 
across  the  search  box  during  the  time  the  satellite  elevation  increases  by  the  height  of  one  camera  field 
of  view  (FOV).  At  the  same  time,  the  camera  elevation  will  move  with  the  nominal  elevation  of  the 
satellite.  The  azimuth  scan  is  repeated,  one  FOV  lower  each  time,  until  the  whole  box  has  been  searched. 

The  real-time  image  processing  during  the  search  basically  consists  of  averaging  the  IR  signal 
received  from  each  pixel  fixed  with  respect  to  the  search  box.  The  processor  tests  for  a  pixel  with  a  signal 
significantly  above  background  after  the  maximum  allowed  integration  time.  The  detailed  nature  of  the 
scan,  whether  step-and-starc  or  continuous  motion,  will  depend  on  the  capabilities  and  limitations  of  the 
available  hardware  and  software.  A  smooth  scan  is  preferred  so  that  time  is  not  lost  waiting  for  the 
telescope  motion  to  settle  before  staring  at  each  new  field.  However,  more  real-time  computational  power 
is  needed  to  integrate  the  signals  from  points  which  appear  to  move  across  the  focal  plane. 

Once  a  probable  satellite  detection  has  been  made  in  one  FOV,  the  telescope  cop*  "oller  is  instructed 
to  stop  scanning  and  new  integrations  are  commenced  on  the  same  field  of  view.  If  the  second  integration 
fails  to  show  a  significant  signal,  a  “false  alarm**  is  declared  and  the  search  resumes  at  the  next  FOV. 
This  practice  of  taking  a  second  look  at  a  suspected  satellite  permits  us  to  lower  the  detection  threshold, 
as  a  much  higher  false-alarm  rate  per  FOV  may  be  tolerated  than  if  only  one  pass  over  the  search  box 
were  permitted. 

If  the  satellite  has  not  been  discovered  by  the  time  the  box  has  been  completely  searched,  a  new 
search  begins,  this  time  with  recomputed  parameters  of  range,  elevation,  search  box  size,  and  integration 
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times.  The  new  search  should  be  more  sensitive  because  the  satellite  is  closer  and  at  a  higher  elevation 
angle  the  second  time. 

4.7  INFRARED  TRACK  AFTER  ACQUISITION 

The  staring  FPAs  are  well  suited  to  tracking  the  satellite  once  it  has  been  acquired.  The  telescope 
would  be  pointed  directly  at  the  detected  satellite,  and  a  tracking  algorithm  determines  the  satellite  image 
centroid  after  each  integration  period.  The  error  signals  thus  obtained  are  used  to  keep  the  system  pointed 
n  the  satellite.  The  angular  uncertainty  in  the  tracking  signals  will  be  roughly  the  angular  size  of  the 
satellite  image  divided  by  the  SNR  of  the  detection.  The  bandwidth  of  the  tracking  signals  will  corre¬ 
spond  to  the  integration  time  in  use.  Further  analysis  of  tracking  is  beyond  the  scope  of  this  report. 
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5.  PERFORMANCE  ANALYSIS 


5.1  RANGE  EQUATION 

The  goal  of  tliis  analysis  is  to  determine  whether  a  given  satellite  can  be  detected  reliably  at  the 
necessary  range  by  a  proposed  system.  Reliable  detection  means  that  the  satellite  signal  S  must  be  larger 
than  the  total  system  noise  by  a  factor  sufficient  to  ensure  a  high  probability  of  detection  Pp  while 
permitting  a  low  rate  of  false  alarms  A  detection  threshold  is  set  at  X  limes  the  noise  level,  where 
X  is  determined  from  K^a  and  the  normal  probability  integral  (assuming  that  the  system  noise  follows 
a  Gaussian  distribution).  The  signal  must  be  above  XN  by  some  amount  YN  to  ensure  a  high  probability 
of  detection  Pp.  Finally,  the  required  SNR  Z^/zv  is  X  plus  Y. 

An  analysis  of  the  required  SNR  is  presented  m  Appendix  C.  There  we  show  that  a  value  for  Z/^fs 
of  6  yields  a  deieciinn  probability  of  95  percent,  while  the  false  alarm  rate  is  12  percent  per  FOV.  As 
discussed  in  Section  4.7,  a  relatively  high  rate  of  false  alarms  is  tolerable  because  the  search  procedure 
allows  us  to  take  a  second  (and  third)  look  at  a  suspected  satellite  detection. 

A  figure  of  merit  often  cited  for  infrared  detectors  is  the  normalized  detct.i .  ’y  r^*  Such  a  quantity 
is  often  useful  in  finding  the  SNR  for  detectors  proposed  for  a  given  application.  However,  the  value  of 
D*  depends  on  the  background  flux,  which  is  responsible  for  a  major  component  of  noise.  In  the  present 
case,  we  prefer  to  calculate  the  photon  background  and  use  a  direct  measure  of  the  other  sources  of  noise, 
such  as  noise-equivalent  input  power  or  zero-flux  readout  noise.  D*  is  discussed  further  in  Appendix  D. 

5.1.1  Signal 

The  total  signal  from  the  satellite,  in  photoelectrons,  may  be  estimated  from 

e  =  -^^^^'A/^^T^(A)Tf(A)G7,(AVA  ,  (8) 


where 

hx  =  spectral  radiant  intensity  (infrared  signature)  of  the  satellite  (W/sier  pm), 
Af  =  collecting  area  of  the  telescope  (m^), 

R  =  range  to  the  satellite  (m), 

he  =  Planck’s  constant  times  the  speed  of  light  =  1.24  eV  pm, 

To  =  average  transmission  of  the  receiving  optics, 
r^(A)  =  transmission  of  the  atmosphere. 
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tf<A)=  transmission  of  receiver  passband  filler, 
quantum  efficiency  of  the  detector, 

A),  A2  =  limiting  wavelengths  of  the  filier-deiecior  combination, 

G  =  photoelectric  gain  of  the  detector  (=  1  for  photovoltaic  devices), 

^INT  -  int-graiion  lime  on  each  pixel  in  ihe  search  box, 
and  all  of  the  important  wavelength  dependencies  are  explicitly  shown. 

5,1.2  Noise 

The  total  system  temporal  noise  N  for  one  integration  period  is  given,  in  electrons,  by 

,  (9) 

where 

ng  =  the  quantum  noise  of  detected  backgrounc'  photons  including  (for  photoconductors)  generation- 
recombination  (G-R)  noise,  and 

n£  =  electronic  noise  from  the  detector,  preamplifiers,  etc. 

The  background  photon  noise  (Poisson  noise)  contains  contributions  from  both  the  sky  radiance  and 
emission  of  the  telescope  optics: 

np^  =  B  +  Bj.  .  (10) 

Quantum  noise.  The  numf«r  of  photoc'-cirons  in  each  unresolved  point  image  produced  by  the 
background  sky  radiation  is  given  by 

B  =  Aj.(or^Q^t,^  .  (11) 

•  02) 

where  the  symbols  are  as  defined  before,  except 
^AX  -  spectral  radiance  of  the  atmosphere,  and 
0)  =  total  solid  angle  of  the  satellite  image. 

The  solid  angle  O)  includes  not  only  the  contributions  of  diffraction,  atmospheric  turbulence,  and 
telescope  aberrations  (adding  up  to  the  optical  image  size  OJ^),  but  also  (Opij ,  the  size  of  the  discrete  detector 
pixels.  The  angular  size  of  the  total  image  results  from  the  convolution  of  the  point-spread  functions  from 
the  individual  contributions,  each  of  which  can  be  characterized  by  a  variance.  The  variance  of  the 
convolved  image  is  the  sum  of  the  individual  variances.  Because  the  variances  vary  as  the  squares  of  the 
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angular  sizes,  they  vary  directly  with  the  solid  angles  of  the  individual  contributions.  By  this  heuristic 
reasoning,  we  may  simply  add  the  solid  angle  contributions 


(13) 

A 

pix 

(0  , 
pix  y2 

(14) 

where 

/  =  the  effective  focal  length  of  the  telescope  and 

^pix  -  the  area  of  one  pixel  in  the  focal  plane. 

The  additional  background  current  caused  by  radiation  of  the  optical  elements 

is  calculated  from 

(15) 

where 

Ljx  =  spectral  radiance  of  a  blackbody  at  the  temperature  of  the  telescope, 
(0  =  solid  angle  subtended  by  a  satellite  image,  as  before,  and 
Csff  =  effective  emittance  of  the  telescope  optics  (Equation  (6)). 

The  total  quantum  noise  power  is  thus  given  by 

~ ^EFF^t]^INT  ’ 


Electronic  nobe.  Sources  of  electronic  noise  include  Johnson  noise  in  the  detector,  shot  noise  from 
detector  dark  current,  noise  from  the  readout  process,  and  preamplifier  noise.  Rather  than  attempting  to 
isolate  and  estimate  each  of  these  sources,  we  have  used  measured  values  of  the  noise-per-pixel  readout, 
obtained  in  conditions  of  minimal  photon  flux  on  the  detectors.  These  noise  values  are  listed  in  Table  !. 
For  the  ground-based  sensors  of  this  study,  we  find  that  the  quantum  noise  from  background  pho»0-.A  is 
the  dominant  term,  in  any  case. 

Many  pixel  readouts  are  summed,  in  general,  to  obtain  a  single  integration  of  satellite  signal.  The 
summation  is  over  individual  frames  which  are  read  during  the  integration  time  and  over  the  number  of 
pixels  which  encompass  the  satellite  image.  The  electronic  noise  is  thus  given  by 


where 


INT 


(18) 


m£  =  number  of  pixels  summed  to  include  the  image  of  the  satellite 


> 
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Np)  =  loiaJ  electronic  noise  (in  electrons)  for  a  single  pixel  readout,  and 
F  =  frame  rate  (sec  *),  assuming  that  ^  2,  3 . 

Total  temporal  r  mse.  If  we  add  together  all  of  these  contributors,  we  find  the  expression  for  the  total 
temporal  noise  is 

^TEMP^  ^EFF^t)'^  R\  W  ' 


Quantization  noise.  In  order  to  perform  the  signal  summations  mentioned  above,  the  analog  pixel 
readouts  are  first  converted  to  digital  form.  The  analog-lo-digital  (A/D)  converter  will  have  enough 
resolution  for  this  task  if  the  least-significant-bil  (LSB)  step  size  A  (here,  in  units  of  electrons)  is  less  than 
approximately  twice  the  RMS  temporal  noise  in  the  analog  signals.  The  RMS  quantization  error  for  a 
single  readout  is 


A 


(21) 


The  quantization  noise  for  the  summed  frames  is  given  by 


—  Ft,.rr 
Q  r  J2  tNT 


(22) 


As  long  as  //v/- »  HF  we  can  treat  the  quantization  noise  as  a  part  of  the  temporal  electronic  noise, 
and  the  total  noise  becomes 


=  N  ^ 
^  ^  TEMP 


+v= 


O^A  ^EFF^' 


F\t 
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Spatial  Soise.  In  addition  to  temporal  noise  and  quantization  noise,  we  must  consider  spatial  noise, 
often  called  fixed-pattern  noise.  Unlike  temporal  noise,  spatial  noise  bears  a  fixed  ratio  to  the  total  signal 
plus  background,  independent  of  integration  time.  Spatial  noise  sets  a  lower  limit  to  the  strength  of  signals 
which  can  be  detected,  even  with  unlimited  integration  time;  it  arises  from  individual  variations  in  the 
responsivity  and  offset  of  each  detector  in  the  array.  In  general,  these  variations  caimot  be  calibrated 
perfectly*^^ 

However,  the  present  application  permits  us  to  minimize  some  of  the  sources  of  spatial  noise.  First, 
the  background  radiance  responsible  for  practically  all  the  photocurrent  should  be  quite  smooth  over  the 
angular  scales  of  interest  here.  Many  frames  of  such  background  can  be  recorded  and  averaged  within 
minutes  of  the  expected  appearance  of  the  satellite.  Furthermore,  this  will  be  done  for  the  same  patch  of 
sky  to  be  searched.  Both  the  magnitude  and  spectral  distribution  of  this  background  radiance  should  be 
unchanged  during  the  acquisition  search.  Thus,  two  of  the  sources  of  spatial  noise  in  the  detectors,  i.e., 
differential  nonlinearities  and  differential  spectral  response,  are  nullified  when  the  average  background 
frame  is  subtracted  from  each  search  frame.  One  source  remains,  which  we  are  unable  to  evaluate  without 
experience  with  the  actual  detector  systems.  This  is  excess  low-frequency,  or  **l/f,”  noise  in  the  detectors. 
This  type  of  spatial  noise  causes  the  irregularities  to  grow  as  the  logarithm  of  the  time  since  the  reference 
frames  were  recorded. 
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To  measure  and  adequately  correct  for  both  the  detector  responsivities  and  the  sky  backgrounds,  the 
camera  signals  must  be  converted  to  digital  form  with  great  accuracy.  Not  only  must  the  quantization 
noise  be  small,  but  irregularities  in  the  signal  thresholds  for  each  digital  output  value  (differential 
nonlineanlies)  in  the  A/D  convener(s)  must  be  kept  to  a  minimum  [10].  At  present,  most  of  the  A/D 
converters  which  can  operate  at  the  speeds  required  for  these  focal  planes  have  a  precision  of  just  12  bits; 
this  probably  accounts  for  reports  that  the  minimum  achievable  spatial  noise  for  HgCdTe  detector  arrays 
is  roughly  0.03  percent  [  1 1  ].  We  wish  to  develop  a  system  in  which  the  the  A/D  errors  are  small  enough 
not  to  introduce  significant  spatial  noise.  Averaged  frames  will  have  to  be  stored  with  more  bits  of 
precision  than  are  provided  by  the  A/D  converters. 

In  infrared  astronomy,  the  procedure  generally  used  to  obtain  accurate  measurements  of  the  flux 
from  a  point  source  is  to  move  a  single  detector  optically  on  and  off  the  source  and  measure  the  difference 
in  signal;  the  background  radiance  is  thus  subtracted  out.  A  staring  sensor  uses  separate  detectors  to  make 
the  on-source  and  off-source  measurements  simultaneously.  If  the  spatial  noise  in  the  array  is  too  great, 
however,  the  search  procedure  could  be  modified  to  scan  the  array  smoothly  across  the  search  box  while 
the  photon  signals  are  integrated  in  computer  memoi^,  as  in  a  time  delay  and  integration  (TDl)  sensor. 

Three  lessons  are  to  be  learned  from  a  consideration  of  spatial  noise.  First,  the  detectors  should  be 
as  uniform  as  possible.  Second,  the  wavelength  passband  of  the  cold  filter  should  be  carefully  chosen  to 
minimize  the  background  radiance,  which  will  be  multiplied  by  the  spatial  noise  factor,  without  unduly 
reducing  the  satellite  signal.  Third,  the  A/D  converters  must  be  accurate  enough  to  obtain  many-frame 
averages  of  high  precision. 

For  the  present  analysis,  we  will  assume  the  spatial  noise  to  be  equal  to  the  residual  temporal  noise 
in  a  one-second  average  of  the  background  radiance.  We  wish  to  emphasize  that  only  experience  with 
the  detector  system  will  show  whether  spatial  noise  can  actually  be  reduced  this  low. 

Clutter  There  can  also  be  real  irregularities  in  the  background,  called  clutter,  which  have  a  similar 
effect  on  the  detection  limits  of  the  system.  In  this  application,  stars  arc  one  form  of  clutter.  From  the 
results  of  the  Infrared  Astronomical  Satellite  (IRAS)  survey  [12],  the  density  of  sources  in  the  12-pm 
band  at  high  galactic  latitude  with  spectral  flux  densities  above  0.5  Jansky  is  approximately  0.6  per  square 
degree  (1  Jansky  =  10*^^  W  m‘^  Hz‘^).  The  number  of  sources  above  a  given  flux  level  decreases  as  the 
-1.0  power  of  the  flux  level.  The  most  sensitive  detector  system  considered  so  far  (HgCdTe)  may  be  able 
to  detect  a  satellite  of  100  W  sier^  at  a  distance  of  4400  km  at  10°  elevation  angle  (in  a  typical  en¬ 
gagement).  This  corresponds  to  a  flux  level  in  the  !2-pm  IRAS  band  of  41  Janskys.  The  corresponding 
density  of  sources  is  0.0073  deg'^  at  high  galactic  latitude,  but  greater  than  0.3  deg‘^  within  a  few  degrees 
of  the  galactic  plane.  A  satellite  position  uncertainty  of  ±  50  km  at  4400  km  leads  to  a  search  box  of  1.7 
deg^ .  The  chance  that  an  IR  star  as  bright  as  the  satellite  will  appear  in  the  box  thus  »-anges  from  -  1 
percent  over  most  of  the  sky  to  ~  50  percent  in  the  galactic  plane.  The  satellite  would  still  be  distinguishable 
by  its  motion  but  the  possible  presence  of  stars  must  be  considered  in  the  design  of  the  search  algorithms. 

Atmospheric  turbulence  already  discussed  in  Section  3.3,  results  in  a  form  of  clutter  called  sky 
noise.  This  is  not  to  be  confused  with  the  Poisson  noise  of  arriving  background  photons.  This  phenom¬ 
enon  arises  from  fluctuations  in  temperature  and  water-vapor  density,  described  by  and  C^.^,  which 
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each  affect  the  local  atmospheric  radiance.  Observations  of  sky  noise  with  individual  detectors  indicate 
that  the  overall  RMS  amplitude  is  about  10*^  of  the  average  sky  radiance.[13]  However,  we  will  only  be 
concerned  with  the  variations  over  angular  scales  of  ~20  to  -100  prad.  Some  of  this  variation  may  also 
be  averaged  out  as  the  satellite  is  tracked  across  the  sky. 

A  very  unpredictable  form  of  clutter  is  from  clouds.  High,  thin  cirrus  clouds,  which  may  not  even 
be  visible  to  the  eye,  can  cause  great  difficulty  in  attempting  to  establish  the  sky  background  level.  In 
addition,  they  attenuate  the  satellite  signal  by  an  unknown  factor.  In  this  report,  all  calculations  assume 
clear  sky  conditions,  i.e.,  no  clouds  present. 

Finally,  other  satellites,  or  even  aircraft,  could  be  confused  with  the  satellite  of  interest.  It  is  assumed 
here  that  as  soon  as  a  candidate  detection  is  made  a  track  will  l>e  established;  it  will  be  rejected  as  clutter 
unless  the  motion  parameters  reasonably  match  those  expected  of  the  satellite. 


5.U  Maximum  Range  of  Detection 


If  we  know  the  system  noise  level  N  and  have  determined  the  required  SNR  Zmin  (=  S/N),  we  can 
use  Equation  (7)  to  find  the  maximum  range  of  detection: 
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In  the  background-limited  case  (ns  »  /*£ ),  we  can  use  Equations  (9),  (17),  and  (6)  to  find 
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It  may  be  observed  from  these  equations  that  the  maximum  detectable  range  is  a  relatively  insen¬ 
sitive  (fourth-root)  function  of  parameters  such  as  telescope  collecting  area,  integration  time,  or  back¬ 
ground  radiance. 


5.2  VALUES  OF  SYSTEM  VARIABLES 

General  system  parameters  used  in  the  performance  calculations  are  summarized  in  Table  2.  Param¬ 
eters  optimized  for  each  type  of  focal  plane  array  are  listed  in  Table  3. 

5.2.1  Signal  and  Noise  Levels 

To  estimate  the  signal  and  noise  levels  for  a  typical  satellite,  we  have  carried  out  the  integrations 
for  Q,  and  Qr,  Equations  (8),  (12),  and  (16).  using  the  characteristics  of  each  type  of  detector  as 
summarized  in  Table  1 .  For  Q,  a  reference  satellite  spectral  radiant  intensity  /sa  was  computed  for  a 
projected-area  emittance  product  of  5  m^,  with  a  temperature  of  300  K.  The  assumed  range  was  1000  km. 
The  background  sky  radiance  for  calculation  of  Qa  is  as  shown  in  Section  3.1  for  clear,  midlatitude 
summer  sky.  Table  3  lists  the  minimum  camera  frame  rates  required  to  avoid  saturation  of  the  detectors. 
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TABLE  2 


Basic  System  Parameters  Used  In  Performance  Calculations 


Telescope  parameters 

Collecting  area 

0.272  m2 

Optical  transmission 

0.74 

Effective  emittance  of  optics 

0.26 

Temperature 

300  K 

Atmospheric  coherence  length  {X  *  0.5  ^im,  zenith) 

5  cm 

LOWTRAN7  parameters 

Atmospheric  model 

Midlatitude,  summer 

Aerosol  model 

Rural 

Visibility  range 

23  km 

Satellite  position  uncertainty  radius 

50  km 

TABLE  3 

Paranteters  for  Different  Fo«^l  Plane  Arrays 


Camera  type 

InSb 

HgCdTe 

SI:Ga 

PtSI 

IrSl 

Effective  focal  length  (m) 

2.5 

2.0 

2.0 

2.5 

1.6 

Relative  aperture  (f/#) 

4.19 

3.35 

3.35 

4.19 

2.68 

Field  of  view  (mrad) 

5.12 

3.20 

3.20 

3.9  X  5.12 

6.40 

FPA  solid  angle  (x  10^  ster) 

26.2 

10.2 

10.2 

20.0 

41.0 

Filter  cut*on  wavelength  (pm) 

3.4 

8.4 

8.4 

3.4 

3.4 

Filter  cut-off  wavelength  (pm) 

4.1 

9.4 

12.0 

4.1 

4.1 

Maximum  pixel  current  (pA) 

47.3 

4681 

3190 

1.36 

3.45 

Minimum  frame  rate  (sec*^) 

10 

1000 

500 

10 

20 

The  detector  signal  currents  resulting  from  this  reference  satellite  are  plotted  as  a  function  of  eleva¬ 
tion  angle  in  Figure  10.  Figure  1 1  shows  the  DC  sky  radiance  currents  vs.  elevation  angle  plus  the  tele¬ 
scope  background  radiance  current.  The  solid-angle  factor  for  these  plots  equals  the  size  of  the  image  plus 
the  size  of  one  pixel.  The  noise  current  density  (in  electrons/Vsec)  caused  by  these  backgrounds  is  shown 
in  Figure  12. 
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Figure  11.  Total  background  currents  vj.  elevation  angle  within  the  solid  angle  subtended  by  an  unresolved 
object. 
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ELEVATION  (d«g) 


Figure  12.  Noise  electron  current  per  triage  area  vs.  elevation  angle  for  various  detectors.  Readout  noise  has 
been  included. 


SX2  Integration  times 

The  integradon  time  is  limited  to  the  amount  of  time  that  the  telescope  can  be  allowed  to  dwell  on 
a  single  FOV  within  the  area  of  sky  to  be  searched.  This  maximum  integration  time  on  any  point  in  the 
search  box  is  the  same  whether  the  telescope  scans  smoothly  across  the  box  or  steps  quickly  from  field 
to  field. 

In  order  to  carry  out  the  satellite  search  at  a  given  elevation,  the  scan  across  the  horizontal  width 
W  of  the  search  box  must  be  completed  iiv  the  time  it  takes  for  the  satellite  to  rise  in  elevation  by  the 
height  Act  of  one  FOV.  The  fraction  of  the  horizontal  scan  which  is  spent  on  any  one  FOV  is  equal  to 
the  width  A0  of  the  FOV  divided  by  IV  (assuming  W  ^  A0).  The  equation  is 

Acr  A0 

- ^  »  (26) 

(da/dt)  W 
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where 

dcjdt  =  the  elevation  component  of  the  angular  velocity  of  the  satellite. 

Tl:e  angular  velocity  of  the  satellite  may  be  calculated  from  its  orbital  parameters  and  its  position 
with  re.spect  to  the  observer.  Assuming  a  circular  orbit,  may  be  found  for  any  position  of  the  satellite. 
The  equations  are  given  in  Appendix  B  for  the  case  of  a  satellite  in  a  circular  orbit  passing  overhead. 
Typical  values  of  as  a  function  of  satellite  position  are  plotted  in  Figure  13.  For  that  figure,  it  is 
assumed  that  the  camera  FOV  is  3.2  mrad  square  and  the  satellite  position  uncenainty  radius  is  50  km. 
For  detection  at  the  ranges  of  greatest  interest,  tfsr  should  be  several  tenths  of  a  second.  Notice  that  the 
background  radiances  will  cause  the  detector  charge  capacities  to  be  filled  in  much  shorter  times.  There¬ 
fore,  to  take  advantage  of  the  possible  integration  times  many  frames  of  pixel  data  will  have  to  be 
summed. 

Notice  that  ifjsjr  depends  directly  on  the  product  AocA(p,  which  is  the  solid  angle  of  the  FOV.  By 
design,  the  solid  angle  of  a  single  pixel  matches  the  optical  image  size  determined  by  turbulent  seeing 
and  diffraction.  Thus,  the  integration  time  is  proportional  to  the  number  of  pixels  (detector  elements)  in 
the  FPA.  This  is  a  general  result  which  holds  for  scanning  systems  as  well  as  purely  staring  systems. 
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Figure  13.  Typical  contours  of  integration  time  vs,  satellite  position.  For  this  plot,  the  relevant  system  parameters 
are  satellite  position  uncertainty  radius  =  50  km,  camera  FOV  =  3.2  x  3.2  mrad. 
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53  MAXIMUM  RANGE  VS.  ELEVATION 


Wc  have  compared  the  performance  of  the  various  FPAs  for  delecting  saielliies.  A  set  of  spreadsheet 
flics  was  set  up  to  find  the  SNR  for  a  given  satellite,  as  a  function  of  satellite  position,  for  each  of  the 
cameras.  The  appropriate  system  parameters,  as  well  as  could  be  determined,  were  used  for  each  different 
focal  plane.  For  simplicity,  all  satellites  were  assumed  to  be  in  circular  orbits  passing  overhead.  Integra¬ 
tion  times  were  taken  to  be  a  function  of  satellite  position  as  described  above.  Two  types  of  plots  have 
been  produced  to  show  the  results.  Figure  14  is  an  X-Y  plot  of  the  positions  of  a  reference  satellite  where 
the  SNR  would  be  6.  The  Y  axis  is  the  observer’s  zenith  direction  and  the  X  axis  is  the  horizontal  distance 
toward  the  rising  satellite.  The  reference  satellite  was  assumed  to  be  at  a  temperature  of  300  K  and  to 
have  a  projecied-area  emiiiance  product  of  5  m^.  This  corresponds  to  a  signature  in  the  LWIR  band  of 
192  W/sier.  The  curves  show  that  the  arrays  with  high  quantum  efficiency  are  capable  of  detecting  such 
a  satellite  at  extremely  long  ranges.  HgCdTe  is  the  best,  with  Si:Ga  a  fair  second;  InSb  is  still  quite  good. 
The  Schotlky-barrier  silicide  detectors  are  relatively  inferior  in  this  comparison.  Nevertheless.  Caniella 
[1]  has  shown  that  a  PtSi  sensor  is  quite  useful  in  the  space  surveillance  application.  His  results  for  the 
maximum  range  of  detection  arc  comparable  to  those  presented  here,  though  they  were  obtained  by 
different  methods. 

5.4  MINIMUM  SATELLITE  SIGNATURE  VS.  POSITION 

The  second  type  of  plot  shows  contours  of  minimum  detectable  signature  as  a  function  of  position  for 
a  given  detector.  Once  again  Zj^i^  =  6.  In  Figures  15-19,  these  threshold  signature  contours  are  presented 
for  the  various  detectors,  in  W/ster  in  the  appropriate  wavelength  band,  3-5  pm  or  8-12  pm.  In  this  way, 
the  results  may  be  used  for  satellites  at  a  wide  range  of  effective  temperatures.  (Sec  Section  2.1  for  the 
relationship  between  signatures  in  the  different  bands.)  Overlaid  on  the  contours  are  circular  orbits  at 
various  altitudes  and  lines  of  constant  elevation  angle.  Such  a  plot  can  be  used  to  answei  a  number  of 
different  questions.  For  example,  given  a  satellite  altitude  and  IR  signature,  the  plots  show  the  minimum 
elevation  angle  at  which  the  satellite  is  detectable  with  each  FPA. 
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Figure  15.  Minimum  detectable  in-band  IR  signature,  as  a  function  of  position,  for  HgCdTe  array  (!28  x  128). 
Signature  units  are  Wfster  in  the  wavelength  band  8-1 2  \im  (LWIR).  Satellites  are  assumed  to  be  passing  overhead. 
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Figure  16.  Same  as  Figure  15.  but  for  Si:Ga  array  (128  x  128). 
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Figure  18.  Same  as  Figure  17.  but  for  IrSi  array  (256  x256). 
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Figure  J9.  Same  as  Figure  17,  but  for  PtSi  array  (160  X  244). 
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6.  DISCUSSION  AND  CONCLUSIONS 


All  of  the  IR  FPA  detectors  we  have  considered  here  show  considerable  promise  for  the  acquisition 
of  satellites  whose  orbits  are  approximately  known.  Predicted  performance  varies  greatly,  from  that  of  the 
PtSi  Schoitky-barrier  array  to  that  of  the  HgCdTe  array.  The  former  could  acquire  medium-sized  satellites 
at  low  altitudes,  especially  at  high  elevation  angles;  the  latter  can  detect  the  smallest  satellites  at  rather 
high  altitudes  or  low  elevation  angles.  All  satellites  are  difficult  to  acquire  at  very  low  elevation  angles. 
This  may  be  seen  in  Figure  14,  from  the  way  all  the  signature  contours  approach  the  origin  ai  low 
elevations.  It  is  interesting  to  observe  that  the  InSb  array  is  ten  times  more  sensitive  than  the  PtSi  array, 
which  also  operates  in  the  MWIR  band.  This  is  to  be  expected,  as  the  InSb  array  offers  roughly  100  times 
the  quantum  efficiency  of  the  Schottky-barrier  detectors. 

Factors  other  than  theoretical  performance  may  also  weigh  in  the  choice  of  FPAs.  For  example,  the 
operating  temperature  of  the  array  affects  the  cost  and  complexity  of  the  cooling  system  required.  We  do 
not  yet  have  sufficient  data  on  the  nonuniformities  or  1/f  noise  of  the  arrays  to  evaluate  spatial  noise, 
though  we  expect  this  to  have  a  significant  effect  on  actual  performance.  The  silicon-based  arrays  are 
expected  to  be  superior  in  this  regard,  if  only  because  the  technology  of  silicon  is  better  developed. 
Finally,  there  is  the  question  of  the  producibility  of  each  array  and  the  percentage  of  usable  pixels  which 
can  be  expected.  This  last  point  is  especially  relevant  to  HgCdTe.  as  many  years  of  development  effort 
by  the  detector  industry  have  not  yet  resulted  in  the  ready  availability  of  large  arrays  of  these  detectors. 

For  the  proposed  experiment,  two  of  the  candidate  arrays  appear  most  promising.  In  the  MWIR 
band,  the  InSb  array  not  only  performs  extremely  well  but  is  readily  available.  Even  better  performance 
can  be  found  in  the  LWIR  band.  Here,  the  Si:Ga  array  looks  very  promising.  Its  theoretical  performance 
does  not  quite  match  HgCdTe,  but  a  working  device  has  already  been  demonstrated  by  Amber  Engineer¬ 
ing;  they  have  achieved  a  high  yield  of  good  pixels  with  good  uniformity.  Although  HgCdTe  could 
operate  at  higher  temperatures  and  show  higher  responsivity,  the  higher  technical  risk  to  achieve  a 
woiidng  system  does  not  seem  to  be  offset  by  these  advantages.  One  complication  of  using  the  LWIR 
band  is  that  the  higher  fluxes  of  background  photons  quickly  saturate  the  charge-holding  capacity  of  the 
pixels.  Thus,  the  arrays  must  be  read  out  at  roughly  1000  frames  per  second.  Digital  circuits  can  be 
designed  to  integrate  the  frames  in  real  lime  so  that  they  can  be  presented  to  the  host  computer  at  a  more 
manageable  rale.  On  the  other  hand,  MWIR  detcaors  do  not  require  frame  rates  in  excess  of  100  sec  ’  to 
avoid  saturation,  and  it  may  be  possible  to  dispense  with  special  frame-integrating  hardware  in  that  case. 

Finally,  altogether  different  arrays  might  be  considered  for  IR  satellite  acquisition;  for  example,  a 
linear  array  might  be  scanned  in  one  direction  to  cover  the  same  search  box.  The  present  analysis  is 
applicable  to  many  of  these  other  cases.  If  the  suggested  detector  has  a  quantum-efficiency  curve  similar 
to  one  of  the  arrays  studied  here,  as  far  as  theoretical  performance  is  concerned  the  main  difference  is 
the  total  number  of  pixels  in  the  array.  The  on-satellite  integration  time  is  proportional  to  this  total 
number  of  pixels  and  the  SNR  varies  as  the  square  root  of  that  number.  For  example,  if  we  substitute 
a  HgCdTe  array  having  480  x  1  pixels  for  one  with  128  x  128  pixels,  the  integration  time  is  reduced  by 
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a  factor  of  34.  the  minimum  detectable  signal  is  increased  by  a  factor  of  5.84,  and  the  range  of  detection 
of  a  given  signature  will  be  reduced  by  the  square  root  of  5.84,  or  2.42.  This  may  still  be  adequate 
performance,  and  the  linear  array  might  offer  offsetting  advantages  over  the  two-dimensional  array. 

The  analyses  contained  in  this  report  show  that  passive  ground-based  acquisition  of  satellites  is 
feasible  with  IR  FPAs  within  the  present  state  of  the  art.  The  remaining  uncertainties,  such  as  the  actual 
magnitude  of  spatial  noise  effects,  can  best  be  resolved  by  building  a  test  system  and  operating  it  in  the 
field.  Such  tests  would  provide  a  base  of  experience  with  different  atmospheric  conditions  while  yielding 
valuable  measurements  of  real  IR  signatures  of  satellites. 
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APPENDIX  A 

EMISSION  OF  OPTICAL  ELEMENTS 


To  estimate  the  thermal  radiation  from  the  optical  elements,  we  set  the  emittance  equal  to  the 
absorpiance  for  the  wavelengths  of  interest.  For  mirrors,  the  absorptance  is  1  minus  the  transmittance.  For 
lenses,  light  is  either  transmitted,  absorbed,  or  reflected.  If  we  assume  that  stray  reflections  from  the  lens 
to  the  detector  come  from  surfaces  at  the  temperature  of  the  lens,  then  the  surface  reflectances  may  be 
added  to  the  absorptance  to  obtain  an  apparent  emittance.  Thus  the  effective  emittance  is  the  complement 
of  the  reflectance  or  transmittance: 

e  =l-r  (A-1) 

(  / 

The  effective  emittance  of  a  series  of  optical  elements  (at  one  temperature)  is  found  from  a  summation 
which  includes  the  transmittances  as  well: 


^EFF  "  ^t  X 


/-t 


^iU^j 


which  reduces  to 


^EFF  "  *  "  fl 


EFF 


(A-2) 


(A-3) 
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APPENDIX  B 

ANGULAR  SPEEDS  OF  SATELLITES 


The  angular  speed  of  the  satellite  depends  on  its  orbit,  its  ra'  e,  and  the  angle  of  the  line  of  sight 
relative  to  the  orbital  velocity  vector.  We  can  simplify  the  problem  by  assuming  that  the  orbit  is  circular 
and  that  the  ground  track  passes  near  the  observer's  position.  Then  for  any  given  cl<*vation  angle  a  and 
range  R,  the  elevation  angular  rate  is  known: 


da 

dt 


1  +“^sina 
R 


dt 


(B-l) 


where 

R  -  slant  range  to  satellite, 

/?£  =  the  radius  of  the  earth,  and 

dyidt  -  the  angular  speed  of  the  orbit  from  the  center  of  the  earth.  In  radians  per  second,  it  is 


dy  _  2k 
dt  ”  5036.9 


where 

//  =  the  altitude  of  the  satellite  orbit. 


{B-2) 
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APPENDIX  C 

REQUIRED  SIGNAL-TO-NOISE  RATIO 


A  search  for  a  weak  signal  in  the  presence  of  noise  must  be  guided  by  two  requirements.  The  first 
is  the  maximum  allowed  rate  (or  probability)  of  false  alarms  Kf^.  The  second  is  the  required  probability 
of  detection  of  a  true  signal  P^.  In  order  to  determine  the  required  SNR  Z^in-  we  also  must  know  how 
many  bins,  or  pixels,  are  to  be  searched  for  the  signal.  Call  that  number  C,  We  assume  that  the  noise  is 
Gaussian  with  standard  deviation  N.  The  signal  S  is  assumed  to  be  small  in  relation  to  the  average  back- 
ground,  so  that  the  bin  containing  the  signal  has  the  same  level  of  noise  as  the  background  bins. 

A  detection  threshold  is  to  be  set  at  a  level  xN  above  the  background.  A  bin  which  exceeds  this  level 
IS  considered  to  contain  a  real  signal.  The  false-alarm  rate  is  computed  as  follows.  The  probability  that 
noise  in  a  single  pixel  will  exceed  xN  is  given  by 


p,w  =  ^[i-vw]  . 

where  is  the  normal  probability  integral 

vw=7^j_y^' • 


For  small  values  of  px  (x  >~3),  it  is  useful  to  use  t  scries  expression: 


,  ,  1  -TVf,  1  1-3  1-3-5 

'  xyjln  L  ^  ^  J 


The  probability  that  no  pixel  will  exceed  xN  is  given  by 


If  pj  «  1 .  we  have 


If  in  addition  pi  «  1/C.  we  have  simply 

K.^CpAx)  . 


The  signal  5  must  exceed  the  detection  threshold  by  an  amount  yN  sufficient  to  ensure  a  high 
probability  of  detection.  The  basic  relationship  is 


In  this  case,  it  is  easiest  to  use  tabulated  values  of  the  nonnal  probability  integral  to  find  a  suitable  y. 
Finally*  the  required  signal  level  is  given  by 

S  =  xN  +  yN  ,  (C-9) 

=  X  +  y  .  (C-10) 

In  the  satellite  acquisition  application,  there  is  a  very  large  number  of  pixels  to  be  searched.  If  the 
range  at  acquisition  is  1 500  km  with  a  position  uncertainty  radius  of  50  km.  then  the  search  box  is  33 
mrad  across.  Typical  image  sizes  of  25  jirad  mean  that  roughly  2,000.000  pixels  need  to  be  searched.  If 
we  set  X  =  5.6,  then  we  find  p\  =  10  *  and  2  percent.  The  false-alarm  rate  is  a  very  steep  function 
of  X  at  this  level.  An  SNR  of  7  allows  v  =  7  -  x  =  1.4,  which  leads  to  Pd  =  92  percent.  Thus,  we  conclude 
that  a  satellite  may  be  acquired  with  high  probability  and  a  low  false-alarm  rate  if  the  SNR  is  7  (-17  dB) 
or  greater.  A  lower  ratio  such  as  6  can  be  used  and  still  give  a  high  Pd*  but  only  if  a  higher  rate  of  false 
positives  can  be  tolerated. 
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APPENDIX  D 

D* 


A  figurc-of-mcrii  often  quoted  for  various  types  of  infrared  detectors  is  the  specific  detectivity, 
written  as  D*.  It  is  defined  as  the  SNR  in  a  1-Hz  bandwidth  per  square  root  of  detector  area.  The  units 
of  D*  are  cm  Hz*^AV.  It  is  related  to  noise-equivalent  power  (NEP)  by 


D*  = 


NEP 


where 

A[)  =  the  area  of  the  detector,  or  pixel  (cm^),  and 
the  bandwidth  (Hz). 

In  turn,  NEP  is  related  to  the  responsivity  R  (A/W),  the  RMS  noise  N  (amperes),  and  the  quantum 
efficiency  T]  by 


NEP  =  ^ 
R 


where  hv  is  the  photon  energy  (eV). 

To  make  use  of  D*,  we  must  also  know  the  conditions  to  which  a  quoted  D*  value  applies:  the  incident 
wavelength  or  spectrum,  the  center  frequency  of  the  temporal  passband,  and  the  total  incident  flux,  or 
solid  angle  in  the  case  of  a  blackbody  spectrum.  Most  frequently,  D*  values  arc  quoted  for  a  295  K 
blackbody  incident  spectrum,  filling  2n  stcradians  in  front  of  the  detector.  Under  such  conditions,  the 
NEP  of  a  good  IR  detector  should  be  dominated  by  the  quantum  noise  of  the  detected  photons.  If  so,  it 
is  said  to  be  a  BLIP  (background-limited  performance)  detector.  The  D*  values  are  useful  in  comparing 
different  detectors  as  long  as  the  conditions  of  the  proposed  application  are  not  very  different  from  those 
for  which  the  values  are  quoted. 

If  the  detector  is  to  be  used  in  a  low-background  situation,  then  the  quantum  noise  will  be  reduced.  For 
example,  reducing  the  FOV  of  the  detector  (increasing  the  fMumber)  will  reduce  the  incident  flux  of 
background  radiance  as  the  sine  of  the  half-angle  of  the  cone  of  incidence.  However,  intrinsic  noise 
sources  may  then  become  dominant.  In  applications  for  which  the  incident  flux  is  low,  or  its  spectrum 
is  not  a  295  K  blackbody,  it  would  be  most  useful  if  D*  values  were  given  for  a  zero  field-of-view,  i.e., 
as  if  a  cold  shroud  completely  covered  the  detector.  It  would  then  be  an  easy  matter  to  obtain  the  system 
noise  floor  of  a  proposed  detector  system  to  compare  with  the  signals  to  be  detected  and  the  quantum 
noise  of  the  expected  photon  background. 
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